Domain-wall and overlap fermions at nonzero quark chemical potential 
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We have recently given a construction of the overlap Dirac operator at nonzero quark chemical 
potential. Here, we introduce a quark chemical potential in the domain- wall fermion formalism and 
show that our earlier result is reproduced if the extent of the fifth dimension is taken to infinity 
and its lattice spacing is taken to zero. We also extend this result to include a bare quark mass, 
consider its continuum limit, and prove a number of properties of the overlap operator at nonzero 
quark chemical potential. In particular, we show that the relation between the anomaly and the 
index of the overlap operator remains valid. 

PACS numbers: 11.15.Ha, 12.38.Gc 



I. INTRODUCTION 

The phase diagram of quantum chromodynamics 
(QCD) as a function of temperature and chemical po- 
tential has been the subject of intense studies over many 
years. It is physically relevant, e.g., for the study of com- 
pact stars, for ultra-relativistic heavy-ion collisions, and 
for the physics of the early universe. The theoretical 
methods that have been employed to investigate the QCD 
phase diagram include model calculations, effective the- 
ories, perturbative studies at high temperature and den- 
sity, lattice simulations, and recently also the AdS/CFT 
correspondence. For an overview of the literature and 
a summary of the current status, we refer the reader to 
Ref. [1]. 

Lattice simulations are the predominant nonperturba- 
tive tool to study QCD from first principles. The case 
of nonzero temperature T can be implemented on the 
lattice without much effort, and therefore the temper- 
ature dependence of many QCD quantities is very well 
understood, see Ref. [2] for a review. This is not true 
for the case of nonzero baryon density or, equivalently, 
quark chemical potential /i. The reason is that at /i 7^ 0, 
the fermion determinant becomes complex so that stan- 
dard importance sampling methods fail. This is an ex- 
ample of the so-called sign problem, see Ref. [3] for a 
thorough discussion of this problem in the context of lat- 
tice QCD. A number of approaches have been invented to 
deal with this problem, such as reweighting along the crit- 
ical line [4], Taylor expansion [5], and analytical contin- 
uation from imaginary fi [6, 7]. Using these approaches, 
attempts have been made to determine the transition line 
Tc{n) and to locate the critical end-point, see Ref. [8] for 
a review. While the results are encouraging, the sign 
problem remains unsolved in principle. 

Most of the lattice results for QCD thermodynam- 
ics have been obtained with staggered fermions, which 
reduce the fermion doubling problem and implement a 
remnant chiral symmetry on the lattice. For simulations 
with less than four staggered flavors, the so-called rooting 
problem has been the subject of some debate, see Ref. [9] 
for a review. We have no intention to enter this debate 
here. Rather, our aim is to investigate how a quark chem- 



ical potential can be implemented in a fermion operator 
that implements an exact chiral symmetry on the lattice. 

In Ref. [10] we showed how this can be done for the 
overlap operator [11, 12]. (For earlier work with a similar 
focus, see Refs. [13, 14].) The resulting operator, Dov(m)i 
contains the sign function of a nonhermitian matrix. We 
also showed that quenched lattice results obtained with 
this operator agree with analytical predictions from non- 
hermitian chiral random matrix theory at /x [15-17], 
see also [18]. Furthermore, the authors of Ref. [19] have 
shown that our construction of Dov(m) yields the correct 
energy density for free fermions. 

In the present paper, we introduce a quark chemical 
potential in the domain- wall fermion formalism [20-23], 
which can be viewed as a particular truncation of the 
overlap operator. In analogy to the well-known result at 
/i = 0, we show that in the limit in which the extent 
of the fifth dimension is taken to infinity and its lattice 
spacing is taken to zero, our earlier result for Dov(a*) is 
reproduced. We also extend this result to include a bare 
quark mass, consider its continuum limit, and prove sev- 
eral properties of I?ov(/^), including the relation between 
the anomaly and the index. 

We should remark that at present, the topic we address 
here may seem to be mainly of theoretical interest since 
lattice simulations with such an operator, especially at 
/i 7^ 0, are numerically much more expensive than those 
with staggered fermions. However, as computers and al- 
gorithms improve, more and more lattice QCD simula- 
tions will be done with overlap and domain- wall fermions. 
As a first step towards such simulations at /i 7^ 0, we have 
already proposed and tested a new iterative method to 
compute the sign function of nonhermitian matrices [24] . 

This paper is organized as follows. In Sec. II we review 
domain- wall fermions at /i = 0, extend the domain- wall 
action to ^ 7^ 0, and consider the limit of this action for 
infinite extent and zero lattice spacing of the fifth dimen- 
sion. After a short side remark in Sec. Ill, we discuss the 
continuum limit of Do^{^) in Sec. IV. In Sec. V we prove 
a number of properties of Z3ov(m) that were stated, but 
not proven, in Ref. [10]. We conclude with a summary 
and outlook in Sec. VI. 
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II. DOMAIN- WALL FERMIONS AT / AND 
RELATION TO OVERLAP OPERATOR 



A. Domain-wall fermions at /i — 

We start with the definition of the Wilson Dirac oper- 
ator in four dimensions, which is given by [25] 

z?w(M) = (4 + M)-i^(7;+ + 7;-) 



i=l 



with 



- 2 (e"^/ + e"''^4" 



(2.1) 



(2.2) 



where M is the Wilson mass, the U G SU(3) arc the lat- 
tice gauge fields with U-„{x) = Ul{x — i'), the "fv are the 
usual Euclidean Dirac matrices, the Wilson parameter r 
has been fixed at r = 1, the 4-d lattice spacing a has been 
set to unity, and for later convenience we have already 
included a quark chemical potential ^,. 

Domain- wall fermions [20-23] are constructed by intro- 
ducing an additional fifth dimension with lattice spacing 
as and extent Ls- The fcrmion fields now have an addi- 
tional index s = 1, . . . , Lg, while the gauge fields remain 
four-dimensional and do not depend on s. At /x = 0, the 
domain-wall fermion action is given by [22, 23] 

= ^(^i.AV'. - - i'sPLi's-l) , (2.3) 

s=l 

where A = a^D^{^ = 0) -|- 1, the chiral projection opera- 
tors Pr and Pl are defined by Pr/l = ^(1 ±75)) and the 
fermion fields satisfy the following boundary conditions 
in the fifth dimension, 

PRi^L^+i = -mPiiVi , ^lV'o = -mPLi^L, . (2.4) 

The quantity m is a bare quark mass parameter. In 
Eq. (2.3), the Wilson mass has to be in the range —2 < 
M < to obtain renormalizablc solutions in the fifth di- 
mension and to avoid the existence of doublers [21, 22]. 

To take the Ls 00 limit, the domain-wall fermion 
action of Eq. (2.3) is supplemented by a pseudo-fermion 
action to cancel divergences due to the heavy fermions 
in the large-Lg limit [23, 26-31]. We use the pseudo- 
fermion action of Ref. [30], which is given by Eq. (2.3) 
with anti-periodic boimdary conditions, i.e.. m, = 1 in 
Eq. (2.4), but in which bosonic fields are used instead 
of the fermionic ones. As the fermion action is only de- 
fined up to a constant normalization factor, we choose to 
multiply the pseudo-fermion action by 1/2. 



B. Domain- wall fermions at yix 7^ 

We define the domain-wall fermion action at /z 7^ to 

be the same action as in Eq. (2.3), except that D„{0) 
is replaced by £'w(a*)- With very minor modifications, 
the arguments of Ref. [21] leading to the bounds on the 
Wilson mass apply to the case of /i ^ as well, and we 
again obtain the requirement —2 < M < 0. 



C. La — > 00 limit of domain-wall fermions at /ix ^ 

For fi = 0, the connection between domain-wall fer- 
mions and the overlap operator in the Ls 00 limit 
has been exhibited in a number of earlier works, e.g., 
Refs. [31-34]. To make the presentation self-contained, 
we now retrace some of the steps taken in these papers, in 
particular Refs. [33, 34], with small modifications suitable 
for our purposes. 

The idea is to introduce successive spinor transforma- 
tions to diagonalize the Dirac operator in the fifth dimen- 
sion and integrate out the fermion fields. To this end, we 
start with the transformation 

= I p^^^ for 1 < a < - 1 , 

i PrXLs + PlXi ioi s = Ls. ^ ^ 

It is straightforward to show that this transformation 
is orthogonal, with Jacobian equal to 1. Substituting 
Eq. (2.5) into Eq. (2.3) with boundary conditions (2.4) 
yields 

Ls _ 

-S5 = J2 i>s{APR - Pl)Xs + MAPr + mPL)xi 

s=2 



+ 2 i's-MPL - Pr)Xs + i>Ls [APl + mPR)xi . 

(2.6) 



s=2 



To simplify the first term in Eq. (2.6), we introduce the 
transformation 



Vi. = Xs{APr - Pl)-^ , 



(2.7) 



which is diagonal in the fifth dimension. The nontrivial 
Jacobian of this transformation can be ignored since it is 
cancelled by the corresponding Jacobian for the pseudo- 
fermions. We also define an operator T by 



T = -{APr - Pl)-\APl - Pr) 



(2.8) 



which is the transfer matrix in the fifth dimension [31] 
and will be discussed in more detail below. Using 
Eqs. (2.7) and (2.8), Eq. (2.6) becomes 



-S5 = '^XsXs + Xi{Pr - mPL)xi 

s=2 

Ls 

- Xs-iTXs - XLsT{Pl - mPfl)xi , (2.9) 



s=2 
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where we have used 



with 



{APn - PlY^APr = Pr , {APn - PlT^Pl = -Pl , 

{APl - Pr)-^APl = Pl , {APl - Pr)-^Pr = -Pr . 

(2.10) 

(See below for comments on the invertibility of APr — Pl 
and APl — Pr-) The structure of Eq. (2.9) suggests to 
transform the Xs for s > 1 according to 



r]s=Xs- Xs-iT , 
with inverse transformation 



(2.11) 



(2.12) 



i=2 



The Xs for s > 1 are transformed according to 

=r?.+T^=+'-^(PL-mPfl)xi . (2.13) 

Both of these transformations have a Jacobian equal to 
1. Inserting Eqs. (2.11), (2.12), and (2.13) into Eq. (2.9) 
leads to 



-S5 = X) + X1-D4X1 



(2.14) 



s=2 



with 



D4 = Pr- mPL - T^^ {Pl - tuPr) . (2.15) 



The 77s and fjs can now be integrated out trivially. Fi- 
nally, we integrate out xi and xi and obtain, together 
with the corresponding contribution of the pseudo-fer- 
mions, 



det 1)4(171) 
det \Di{l) 



det Deff 



with an effective 4-d operator DeS given by 



D^s = (1 + m) + (1 - m)-^" 



(2.16) 



(2.17) 



We now take a closer look at the transfer matrix of 
Eq. (2.8), which can be rewritten as 



T = (1 + a^H^PR)-\l - a^H^PL) 



(2.18) 



with = 75-Dw For /x = 0, is 75-hermitian, i.e., it 
satisfies = 75-Dw 75) and thus is hermitian. From 
this it follows that T is also hermitian. The transfer 
matrix can be related to a 4-d Hamiltonian Ht by writing 
it in the form 



T 



1 - a^Ht 
1 + a^Ht 



(2.19) 



Ht = {2 + a^H^j^y^H^ = H^{2 + a^l^H^)'^ . 

(2.20) 

For /i = 0, Ht is hermitian. 

Up to this point, everything went through as in 
Refs. [33, 34]. We now move; on to the case of /i 7^ 0, in 
which Dw ceases to be 75-hermitian. Therefore, neither 
i?w nor Ht nor T are hermitian. To obtain the is — > 00 
limit of Eq. (2.17), we consider the matrix function 



f{Ht) 



1 - T^- 



(2.21) 



where T is given by Eq. (2.19) with a nonhcrmitian ma- 
trix Ht- A function / of an arbitrary complex matrix C 
can be defined by [35] 

fiC) = ^J^dzf{z)[z-C]-\ (2.22) 

where the integral is defined component-wise and F is a 
collection of closed contours in C such that / is analytic 
inside and on F and such that F encloses the spectrum 
of C. We are therefore interested in 



a = 



1 



where 



t = 



1 - z 
l + z 



(2.23) 



(2.24) 



and z G C is an eigenvalue of a^Ht- If \t\ < 1 {\t\ > 1), 
(7 — > 1 (fj — > —1) as Lg — > 00, i.e., we can write a — > 
e(l — |tp), where e denotes the sign function. For z G R, 
\t\ < 1 {\t\ > 1) if 2 > (2; < 0), and hence a = e{z). Let 
us now consider z gC with a; = Re 2;, for which 



Ax 



|1 + ^P ■ 

Prom this expression we obtain immediately 
lim a = s{x) = e(Re z) =: s{z) , 

Ls — *oo 



(2.25) 



(2.26) 



where the last equality defines the sign function of a com- 
plex number. This can also be written as 



e{z) 



(2.27) 



with the branch cut of the square root along the negative 
real axis. Wc thus obtain 

lim D^s{tx) = {l + m) + {l-m)j5£{Ht{fi)), (2.28) 

Lg—^oo 

where the sign function £{C) of a nonhermitian ma- 
trix C is defined formally by Eq. (2.22) in combination 
with Eq. (2.26) or (2.27). A simpler form for s{C) can 
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be obtained if C can be diagonalized, i.e., if it can be 
written in the form UMJ^^ with U G Gl(iV,C) and 
A — diag(Ai, . . . , \n), where is the dimension of C. 
It then follows from Eq. (2.22) that 



f{C) = Uf{A)U-' 



(2.29) 



with 



/(A) - diag(/(A.)) (2.30) 

so that the matrix sign function can be defined by [36] 

e{C) = Us{ReA)U-^ . (2.31) 

Note that if C cannot be diagonalized, one can use the 
Jordan canonical form instead, see Ref. [24] for details. 

In the course of the derivation, we have assumed (i) 
that the operators APr — Pl and APl — Pr arc invcrt- 
ible and (ii) that the elements of the diagonal matrix A 
are not zero or purely imaginary so that the sign function 
is well-defined. For any of these assumptions to be vio- 
lated, the gauge-field would have to be fine-tuned. This 
happens on a gauge field set of measure zero and can 
therefore be ignored in practice. The same remark ap- 
plies to the possibility that Ht {ji) might not be diagonal- 
izable. 

D. as — > limit 

To recover the standard overlap operator, it remains 
to take the limit 05 ^ in Eq. (2.28), and therefore in 
Eq. (2.20), which yields 

-Dov(At) = lim„ lim Des{^Ji) 

= (1 + m) + (1 - m)75 e{H^{ti)) . (2.32) 

For m = 0, Eq. (2.32) together with Eq. (2.31) agrees 
with our earlier result [10]. For m 7^ 0, we see that the 
quark mass is included in the same way as for = [31]. 

III. SIDE REMARK 

Note that for ji = 0, there is an expression for the 
overlap operator that is equivalent to Eq. (2.32), i.e., 



£)ov = {l+m) + {l-m) 



(3.1) 



One could be tempted to use this expression for /z 7^ as 
well, with £>w replaced by D„{fi). However, the resulting 
operator is not equivalent to Eq. (2.32) for /i 7^ due 
to the lack of 75-hermiticity of D^{(ji). In particular, it 
does not satisfy the Ginsparg- Wilson condition and has 
no exact zero modes at finite lattice spacing. Thus, the 
expression (3.1) is not suitable for an extension to /x 7^ 0. 



IV. CONTINUUM LIMIT OF r>ov(M) 

In this (and only this) section we reintroduce the 4-d 

lattice spacing a that was set to unity earlier and write 
the Wilson Dirac operator, with the argument fi sup- 
pressed, in the form [37] 



= - [-{l + s)+ab^ 



(4.1) 



where 1 + s = —Ma with jsj < 1 and is the massless 
Wilson operator, = D^{M = 0). We therefore have 

for H„ = 75 -Dw 



{aH^f = {l + sf -a{l + s){b^ 
= (l + s)2+0(a2). 



■ 75^w75) + 0{a^) 
(4.2) 



where the last step follows from the fact that an- 
ticommutes with 75 up to terms of order a. Using 
Eqs. (2.22) and (2.27), the matrix sign function can be 
written as 



and we find 

75e(ifw) = -1 + 



(4.3) 



+ 0(a2) . (4.4) 



1 + s 

Equation (2.32) thus becomes 

Do■v{^J) = ^[(1 +ma) + ma)75£(i?w)] 



= 2m + 



1 + s 



£»w(m) + 0{a) . 



(4.5) 



Since for a — > the Wilson Dirac operator becomes the 
continuum Dirac operator, this is also true for _Dov(m)i 
up to a normalization factor. There was no need in the 
above derivation to use the 75-hermiticity of b^, which 
is lacking for fj, ^ 0. The only input required was the fact 
that Dw anticommutes with 75 in the continuum limit, 
which holds for 7^ as well. 



V. PROPERTIES OF DqvCm) 

In Ref. [10] we stated, but did not prove, a number of 
properties of Z)ov(m) at m = 0. We supply the missing 
proofs here, assuming m = throughout this section. 

Property 1: For n ^ 0, £'ov(m) is no longer 75- 
hermitian but satisfies 

75£>ov(m)75 = Dli-tx) (5.1) 
instead. From Eq. (2.32) with rn = we obtain 

75-Dov(/")75 = 1 + £(i/w(/U))75 , 

Z)t^(-M) = l + £^(Fw(-M))75, (5.2) 
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from which Eq. (5.1) follows because of 

= £(75i5w(M)) = e{H^{ix)) . (5.3) 

In the second step, wc used the well-known fact that D„ 
satisfies Eq. (5.1), which follows from Eq. (2.1) and from 
75T^75 = (T'jF)t. In the third step, we used the fact that 
the sign function satisfies e(C^) = £^(C') for any matrix 
C, which follows from Eq. (2.22) with s*{z) = s{z*), or 
from Eq. (2.31). Note that the proof of Property 1 also 
goes through for m ^ 0. 

Property 2: L'ov(m) satisfies a Ginsparg- Wilson rela- 
tion [38] of the form 



{£),75} =£'75-0 



(5.4) 



Setting m = in Eq. (2.32) and using the shorthand e 
for e{Hy,{n)), we have 

{Dovifi), 75} = 275 + e +75 £75 , 
-Dov(/")75£'ov(m) = (1 + 75 £)75(1 + 75 £) 

= 75 + 75£^+e+75£75 • (5.5) 

Equation (5.4) now follows immediately from £^(C) = 1 
for any matrix C, which in turn follows from Eq. (2.22) 
with e'^{z) = 1, or from Eq. (2.31). 

Property 3: All eigenvalues of £'ov(a*) that are not 
equal to or 2 come in pairs A (with eigenvector ip) and 
A/(A — 1) (with eigenvector 75V'). We first note, with all 
arguments suppressed, that 



implies 



DovV' = (1 + 75 e)V' = 



eijj = {X- 1)75^ , 



(5.6) 



(5.7) 



from which wc conclude that 



£'ov(75V') = 75V' + 75 £ 75V' = 75V' + 75 £ 



£ V' 

A- 1 



A 



A - 1 



(75V') , 



(5.8) 



where we have again used £-^ = 1. Note that for /x = 0, all 
eigenvalues lie on the circle with radius 1 and center at 
1, in which case A and A/(A — 1) are complex conjugates 
of each other. 

Property 4: The mapping A ^ z = 2A/(2 — A) projects 
the pair A and A/(A— 1) to a complex conjugate pair ±z. 
This follows from elementary algebra. 

Property 5 makes statements about the eigenvectors 
of Dov{ii) corresponding to eigenvalue or 2. In these 
cases. Property 3 implies that tp and j^t/j are degenerate 
eigenvectors of Do^{ji). This means that 75 commutes 
with Dov{ij) in the corresponding degenerate subspace 
and can thus be diagonalized in this subspace. Because 
of 75 = 1 the eigenvalues of 75 are ±1, i.e., the eigen- 
vectors of -Dov(a*) corresponding to A = or 2 can be 



arranged to have definite chirality. In the following we 
denote by the number of eigenvectors corresponding 
to A = or 2 with (75) = ±1. Consider now the opera- 
tor B = Dov + JsDov'Jb, where the argument /x has been 
suppressed. It is easily shown that if is an eigenvector 
of with eigenvalue A, then tpx and J5^p\ are degen- 
erate eigenvectors of B with eigenvalue A^/(A — 1). For 
A 7^ 0, 2, we now construct the vectors ^p'f = ip\±j5ip\- 
According to Property 3, ipx and "f5ipx are linearly inde- 
pendent in this case, and therefore the two vectors ipx 
are nonzero and linearly independent. Moreover, they are 
also eigenvectors of 75 with eigenvalue ±1, respectively. 
We now consider _B in a basis consisting of the and 
of the eigenvectors of Z^ov corresponding to A = and 2 
with definite chirality. Since in this basis the operators 
75 and B are simultaneously diagonal with eigenvalues 
A^'*'^' and Ap\ respectively, we have 

tr(75i?)=53Al'^^Up^ 

i 

1 A^ 
= 2 5 ^'^^ ~ '^'^ aT^ ^ ~ ^° ^ ^""^ ~ ""^^^ 

(5.9) 



4(nJ - n^) , 



where in the second line di is the (accidental) degeneracy 
of the eigenvalue Ai 7^ 0, 2 of Dov and the factor of 1/2 
in front of the sum removes a double counting of eigen- 
values. An analogous argument holds for Doy = 2 — Dov, 
which also satisfies the Ginsparg- Wilson relation and for 
which the roles of A = and A = 2 are interchanged, 
leading to 



tr(75S) = 4(n+ - tIq ) , 



(5.10) 



where B = Dov + IsDav^n- From tr(75B) = -tr(75B) 
we conclude that 



nA - n. 



= -("2" - "2 ) 



(5.11) 



as stated in Ref. [10]. From tr(75B) = 2tr(75Dov) we 
also conclude that 



■tr(75Dov) = 2(n(j 



2indcx(i?ov). (5.12) 



The relation (5.12) between the anomaly and the index 
of £>ov was already proven for /x = in Refs. [29, 39, 
40] . Our simple derivation shows that it remains valid at 
/i ^ 0. (The method introduced in Ref. [40] also works 
at /u ^ without modifications.) Eq. (5.12) was used in 
Ref. [10] to explain an observed shift in the number of 
zero modes of £'ov(m) as a function of ji. 

Property 6 concerns the normality of Dovifj)- Prom 
Eqs. (5.1) and (5.4), one easily shows that 



£)ov(m)£>Jv(-m) = DM + Dl{-n) 
= Dl{-^x)Dov{fJ,) . 



(5.13) 



This means that for /i = 0, Dov is a normal operator, 
whereas for ^ 0, we cannot conclude anything from 
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Eq. (5.13) about the normality of Dov{^^)- This suggests 
that for fi 0, D^y^ii) is not a normal operator (at least 
generically) . This expectation is confirmed numerically. 
It is interesting to note that the operator (5.13) is equal 
to the operator B we defined in the proof of Property 5. 

VI. SUMMARY 

We have extended the domain-wall formalism to non- 
zero quark chemical potential and have shown that in 
the limit in which Lg ^ oo and 05 ^ we obtain an 
expression for the overlap Dirac operator that is identical 
to our earlier result [10]. We have also included a bare 
quark mass, considered the continuum limit, and proven 
a number of analytical properties of this operator. 

In actual lattice simulations, the use of Eq. (2.32) will 
be hindered by two problems. The first is the infamous 
sign problem that plagues lattice QCD at 7^ 0. We 
have nothing new to say about this problem here. The 
second problem is that not much is known about efiicient 
numerical computations of the sign function of a nonher- 



mitian matrix. As remarked earlier, we have started to 
address the second problem in Rcf. [24]. 

Work is in progress in several directions. First, we 
will continue our algorithmic developments to compute 
the sign function of nonhermitian matrices, with partic- 
ular emphasis on novel deflation schemes. Second, we are 
currently testing the predictions of nonhermitian random 
matrix theory also for the unquenched theory, which can 
be done by reweighting on the small lattices that we have 
studied so far. We are also studying the average phase 
factor of the fermion determinant, for which some an- 
alytical predictions exist in the epsilon-regime of QCD 
[41, 42]. 
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